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ABSTRACT 

Middle  ultraviolet  spectra  of  the  atmospheric  airglow  were  obtained  from  a  March 
1992  rocket  flight  of  the  NPS  MUSTANG  instrument.  These  spectra  are  analyzed  from 
1900  A  to  3100  A,  over  an  altitude  range  of  100  km  to  320  km.  The  data  are  modeled 
with  computer  generated  synthetic  spectra  for  the  following  emissions:  N2  Vegard  Kaplan 
(VK);  N2  Lyman-Birge-Hopfield  (LBH);  and  NO  Gamma  (y),  Delta  (5),  and  Epsilon  (e) 
bands.  A  best  fit  procedure  was  developed.  The  resulting  synthetic  spectra  agree  well 
with  obtained  airglow  data.  Confirmation  was  made  of  the  theoretical  self  absorption 
versus  non-self  absorption  processes  of  the  NO  (0,0),  (1,0),  (2,0)  y  resonance  band 
emissions.  NO  self  absorption  is  a  necessary  inclusion  of  any  atmospheric  nitric  oxide 
analysis  stratagem.  Profiles  of  temperature  versus  altitude  and  NO  column  density  versus 
altitude  for  the  rocket  flight  are  estimated. 
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I.        INTRODUCTION 

The  electromagnetic  properties  of  the  earth's  atmosphere  are  important  to  many  areas  of 
science  such  as  weather  forecasting,  atmospheric  instruments,  satellite  electromagnetic 
surveillance,  and  communication  operations.  The  upper  atmospheric  region  known  as  the 
ionosphere  contains  electron  densities  which  are  sufficiently  high  to  form  a  plasma  medium. 
The  plasma-like  properties  of  the  ionosphere  are  of  specific  interest  to  the  military.  They 
strongly  affect  operations  of  systems  such  as  high  frequency  (HF)  radio  communications,  over- 
the-horizon  radar  (OTHR),  ballistic  missile  early  warning  (e.g.,  SDI),  and  the  Ground  Wave 
Emergency  Network  (GWEN).  For  effective  prediction  of  the  plasma-like  properties  of  the 
ionosphere  an  electron  density  profile  is  required. 

An  electron  density  profile  is  a  plot  of  electron  density  versus  altitude  for  a  given  time  and 
location.  In  1986,  the  Joint  Chiefs  of  Staff  prioritized  the  determination  of  electron  density  as 
fifth  on  a  list  of  43  critical  geophysical  parameters  (The  Joint  Chiefs  of  Staff,  1986).  The 
electron  density  varies  with  time  (e.g.,  day,  night,  solar  cycle),  altitude,  and  geographic 
location.  Presently,  electron  density  profiles  are  determined  using  approximately  20  ground- 
based  ionosonde  stations  worldwide.  This  can  be  very  effective  locally,  but  is  inherently 
imprecise  for  a  global  forecasting  model. 

An  alternative  to  ionosonde  ground  stations  is  to  deduce  density  profiles  of  the  neutral 
atoms  and  molecules  in  the  ionosphere  from  naturally  occuring  airglow  emissions.  From  these 
neutral  density  profiles  and  knowledge  of  atmospheric  photo-chemistry  one  can  infer  electron 
and  ion  density  profiles.    The  inference  of  electron  densities  from  emission  observations  is 


performed  by  relating  measurements  of  natural  ultraviolet  emissions  to  the  neutral  species' 
densities. 

An  operational  ultraviolet  spectrograph  that  is  relatively  small  and  low-powered  can  easily  be 
incorporated  into  an  orbiting  tactical  communication  platform  or  an  optical  surveillance  satellite 
configuration.  It  would  achieve  continuous,  passive  measurements  of  the  ultraviolet  radiation 
emitted  by  the  ionospheric  airglow.  A  low  earth  orbit  (LEO)  based  satellite  would  provide  a 
globally  complete  remote  sensing  technique  of  predicting  an  electron  density  model  of  the  Earth's 
ionosphere. 

The  Naval  Research  Laboratory  (NRL)  and  The  Ionospheric  Physics  Group  at  the  Naval 
Postgraduate  School  (NPS)  have  a  joint  spectrographic  project.  The  NRL  project  is  a  High 
Resolution  Airglow  and  Aurora  Spectrograph  (HIRAAS)  with  an  observation  range  of  500  A  to 
1500  A.  The  addition  of  the  HIRAAS/MUSTANG  experiment  represents  a  significant  progress  of 
the  global  determination  of  electron  density  profiles  beyond  ground-based  ionosondes.  A  sounding 
rocket-launched  experiment  was  first  successfully  tested  in  March,  1990  from  White  Sands  Missile 
Range,  New  Mexico.  The  data  for  this  thesis  is  from  the  second  successful  launch  which  was 
conducted  in  March,  1992.  A  third  rocket  experiment  will  be  conducted  in  March,  1994,  from  the 
Poker  Flat  Missile  Range  in  Alaska.  The  relative  high  latitude  of  the  Poker  Flat  Range  will  enable 
the  observation  of  the  effects  of  auroral  electron  precipitation  on  ionospheric  photochemistry. 

These  sub-orbital  rocket  experiments  will  be  followed  by  a  remote  sensing  satellite  project.  The 
successor  to  the  MUSTANG  instrument  is  the  Ionospheric  Spectroscopy  and  Atmospheric 
Chemistry  (ISAAC)  instrument  which  is  already  constructed.  The  NPS  Ionospheric  Physics 
Group  is  presently  calibrating  ISAAC.  It  will  be  launched  on  the  Advanced  Research  and  Global 
Observation  Satellite  (ARGOS)  in  late  1995. 


A.  Thesis  Objectives 

This  thesis  focuses  on  the  analysis  of  spectral  data  obtained  from  the  March,  1992  rocket 
flight.  The  analysis  extends  methods  used  previously  on  subsets  of  the  data  to  the  observed  range, 
1900  A  to  3100  A  (see  e.g.,  Cleary,  1986;  Clayton,  1990;  Anderson,  1990;  Mack,  1991;  Walden, 
1991).  All  known  emissions  in  this  range  are  used  in  the  fitting  routines  except  OI  2972  A,  Oil 
2470  A,  OIII  2853  A,  and  Nil  2143  A  atomic  lines.  The  result  of  the  fitting  algorithm  is  a  set  of 
relative  scale  factors  as  a  function  of  altitude.  The  relative  scale  factors  are  converted  into 
absolute  terms  and  plotted  against  altitude.  These  absolute  plots  of  intensity  and  column  density 
can  be  used  to  infer  an  electron  density  profile  (see  e.g.,  McCoy,  1985;  Meier,  1991). 

Additional  goals  of  this  thesis  are  to  confirm  the  theoretical  self  absorption  mechanisms  of  the 
NO  (0,0),  (1,0),  (2,0)  y  band  emissions,  determine  the  NO  column  density  profile  and  estimate  a 
temperature  profile  for  the  observed  altitude  range. 

B.  Thesis  Outline 

This  thesis  is  separated  into  seven  chapters.  Chapter  II  gives  general  background  information 
on  the  descriptions  and  designations  of  the  atmosphere  and  airglow  emissions.  A  basic  explanation 
of  molecular  spectra  is  also  outlined  in  this  section.  Chapter  III  addresses  the  specific  dayglow 
emission  characteristics  in  the  wavelength  range  of  interest.  Chapter  IV  presents  an  overview  of 
the  experiment.  Chapter  V  gives  a  full  breakdown  of  the  algorithm  techniques  to  generate  a 
synthetic  model  to  fit  the  data.  Chapter  VI  explains  the  results  of  the  study.  Chapter  VII 
concludes  the  document  with  a  summary  of  findings  and  suggestions  for  future  study.  The 
appendices  contain  the  plots  of  the  final  fits  to  the  data  as  well  as  an  estimated  NO  column  density 
profile  and  atmospheric  temperature  profile. 


II.       BACKGROUND 


A.    The  Atmosphere 

1.     General  Description 

The  Earth's  atmosphere  is  designated  by  distinct,  horizontal  layers.  These  layers  are 
determined  by  one  of  two  primary  methods.  One  method  designates  the  layers  based  upon 
temperature  characteristics  and  the  other  uses  constituent  density  to  classify  the  atmospheric 
boundaries.  Both  schemes  are  shown  in  Figure  II- 1  .  The  figure  shows  two  temperature 
profiles  corresponding  to  low  and  high  solar  activity. 
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Figure  II- 1    Atmospheric  temperature  profile.    Also  shown  are  atmospheric  layers  designated 

by  temperature  compared  with  corresponding  density  layers. 


This  method  of  classification  designates  the  troposphere  as  the  region  closest  to  the 
Earth's  surface.  In  the  troposphere  temperature  decreases  with  altitude,  on  average  from  about 
300  K  at  the  surface  to  220  K  at  10  km.  The  point  at  which  temperature  begins  to  rise  is  called 


the  tropopause  and  it  determines  the  start  of  the  next  region,  the  stratosphere.  The  stratosphere 
is  a  relatively  warm  layer  due  to  absorption  of  solar  radiation  by  ozone  (primarily  2000-3000 
A).  The  top  of  the  stratosphere  is  marked  by  the  stratopause.  The  next  layer  is  labeled  the 
mesosphere.  In  the  mesosphere,  the  temperature  decreases  with  altitude  until  a  low  of  about 
100  K  occurs  at  approximately  85  km.  This  is  the  mesopause,  which  marks  the  beginning  of 
the  next  layer,  the  thermosphere.  Temperature  rises  steadily  in  the  thermosphere  with 
increasing  altitude  until  reaching  an  equilibrium  temperature  between  1000  K  and  2000  K, 
depending  on  solar  activity. 

Studying  the  atmosphere  in  terms  of  density  or  chemical  composition  gives  rise  to  a 
different  series  of  atmospheric  layers.  The  lowest  80  kilometer  layer  of  the  atmosphere  is 
called  the  homosphere.  Vertical  mixing  of  this  region  provides  a  uniform  relative  density  ratio 
between  various  molecular  species.  The  primary  components  are  molecular  nitrogen  and 
oxygen  and  atomic  argon  in  the  well  known  ratio  of  .78:  .21:  .009.  The  remaining  0.001  is 
primarily  helium,  hydrogen,  and  carbon  dioxide.  The  layer  from  80  km  to  100  km  is 
designated  the  barosphere  and  it  is  a  transition  region  above  which  very  little  vertical  mixing 
occurs.  Atop  the  barosphere  is  the  heterosphere  from  1 00  km  to  approximately  1 000  km.  In 
this  area  of  the  upper  atmosphere,  the  relative  densities  vary  with  altitude  due  to  differing 
molecular  weights.  The  density  of  the  atmosphere  above  the  heterosphere  is  so  low  that 
individual  molecule  trajectories  are  essentially  parabolic  orbits.  This  region  is  labeled  the 
exosphere. 

2.     The  Ionosphere 

The  data  analysis  in  this  work  pertains  to  the  lower  heterosphere.  This  corresponds  to 
the  ionosphere.  Specifically,  the  first  200  km  of  the  thermosphere  where  the  convective  mixing 


is  minimal  and  there  exists  a  positive  temperature  gradient.  The  ionosphere  is  so  designated 
due  to  the  ionization  of  neutral  atmospheric  constituents  into  ions  and  electrons.  The 
ionosphere  is  divided  into  four  layers  due  to  the  presence  of  discrete  electron  density  regions 
and  is  depicted  in  Figure  11-2  .  The  layers  are  labeled  D,  E,  Fi  and  F2.  The  actual  height  and 
thickness  of  each  layer  varies  with  solar  cycle  activity.  The  D  region  is  the  lowest  layer  and 
typically  extends  from  about  50  km  to  90  km.  The  E  region  is  from  90  km  to  130  km.  Fi 
ordinarily  exists  from  130  km  to  160  km  during  daytime  periods  and  the  F2  region  lies  above 
160  km. 
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Figure  JI-2  Atmospheric  electron  density  profile. 

The  primary  mechanism  for  the  development  of  the  ionosphere  during  the  daylight 
hours  is  photo-ionization.  The  general  chemical  equation  below  for  a  single,  neutral  molecule 
M  describes  this  process. 


M  +  hu-^  M+  +e~ 


The  hv  term  represents  solar  photons  with  sufficient  energy  for  ionization.  The  ionization  rate 
Qx(z)  is  proportional  to:  1)  the  flux  of  photons  as  a  function  of  wavelength  7x.(z)  2)  the 
density  of  neutral  molecules  w(z)  exposed  to  the  solar  flux;;  3)  the  efficiency  with  which  these 
photons  cause  ionization.  This  efficiency  is  termed  the  ionization  cross-section  a^X).  It  is 
wavelength  dependent  and  depends  on  the  ionization  potential  of  the  molecule  M.  Thus, 

Qx(z)  =  n(z)*Ii(z)*ai(X) 
The  photon  flux  decreases  in  intensity  with  decreasing  altitude  and  the  neutral  species 
concentration  decreases  with  increasing  altitude.  The  product  of  these  two  quantities  produces 
a  profile  with  a  discrete  layer.  This  profile  is  referred  to  as  a  Chapman  ionization  profile.  This 
is  depicted  in  Figure  II-3.  In  region  (A)  of  this  figure  the  species  concentration  is  not  high 
enough  to  cause  significant  photon  absorption.  Maximum  ionization  occurs  at  a  distinct 
atmospheric  area  (B).  The  ionization  rate  falls  off  rapidly  at  lower  altitudes  because  very  few 
ionizing  photons  are  able  to  penetrate  through  to  the  (C)  section. 
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Chapman  (1931)  was  the  first  to  provide  a  quantitative  explanation  for  the  formation 
of  an  atmospheric  ionized  layer.  This  ionization  process  is  effective  enough  to  maintain 
sufficient  quantities  of  ionized  species  to  maintain  an  ionospheric  plasma.  The  collective 
behavior  of  the  ionospheric  plasma  can  affect  electromagnetic  waves  propagating  through  the 
ionosphere  and  is  of  special  interest  to  the  military. 
3.     Atmospheric  Airglow 

The  MUSTANG  instrument  was  designed  to  measure  emissions  of  the  Earth's 
airglow.  Airglow  is  a  naturally  occurring  radiation  from  the  Earth's  atmosphere.  Dayglow,  or 
daytime  airglow  is  primarily  driven  by  photoexcitation.  Also  contributing  to  dayglow  are 
emissions  from  photoelectron  impacts  and  photochemical  reactions.  The  intensity  of  emissions 
resulting  from  photochemical  reactions  is  negligible  during  daylight  hours  compared  with  the 
intensity  from  photoexcitation.  Photochemical  reactions  are  significant  at  night  since  there  are 
obviously  few  direct  solar  photons.  Photochemical  reactions  are  not  specifically  addressed  in 
this  analysis  because  the  data  was  obtained  during  the  daylight  period. 

Photoexcitation  is  the  excitation  of  an  atom  as  a  result  of  the  absorption  of  a  photon. 
The  atom  drops  back  to  a  lower  energy  state  by  emitting  a  photon  with  energy  equal  to  the 
transition.  This  process  will  be  further  described  in  the  next  section  on  molecular  emissions. 

B.    Molecular  Spectra 

Molecular  emissions  may  occur  as  the  result  of  electronic  state  changes,  but  they  may  also 
occur  from  purely  vibrational  and  rotational  changes  of  state.  Figure  II-4  illustrates  an 
electronic  transition.  Each  of  these  iypes  of  transitions  has  typical  energy  gaps,  with  electronic 
being  on  the  order  of  10  eV,  vibrational  of .  1  eV,  and  rotational  of  .001  eV. 


As  seen  in  the  figure,  for  any  given  electronic  energy  level  there  are  many  vibrational 
energy  levels  labeled  v*  for  the  excited  electronic  state  or  v"  for  the  lower  or  ground  electronic 
state.  Similarly,  for  a  given  change  in  vibrational  level  there  are  many  possible  initial  and  final 
rotational  levels.  As  a  result,  molecular  emissions  from  one  vibrational  to  another  are 
composed  of  many  closely  spaced  lines.  These  emissions  are  called  molecular  bands. 

The  principal  photon-molecule  interactions  are  shown  in  Figure  H-5.  For  elastic 
scattering,  an  energy  transfer  between  the  photon  and  the  molecule  does  not  take  place,  even 
though  the  photon  may  change  direction.  With  inelastic  scattering,  the  incident  photon  not  only 
changes  direction,  but  imparts  energy  to  the  molecule.  In  this  case,  the  molecule  is  excited  to  a 
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Figure  11-4  Molecular  energy  versus  internuclear  distance. 


higher  energy  state.  If  the  energy  difference  between  states  is  equal  to  a  level  AE,  then  the 
deflected  photon  will  propagate  with  an  energy  AE  less  than  the  initial  photon's  incident 
energy.  Another  photon  with  this  energy  AE  will  be  emitted  when  the  molecule  returns  to  the 
ground  state.  Resonance  scattering  is  initiated  when  the  incident  photon  is  absorbed,  raising 
the  molecule  to  an  excited  state.  This  is  followed  by  a  re-emission  when  a  photon  with 
wavelength  equal  to  the  incident  photon  is  emitted.  Fluorescence  scattering  occurs  when  the 
excited  molecule  undergoes  several  downward  L-isitions  in  returning  to  the  ground  state. 
Each  transition  generates  a  separate  photon.  An  excited  molecule  can  also  collide  with  another 
particle,  thus  returning  to  the  ground  state  without  emitting  a  photon.  This  is  called  collisional 
de-excitation  or  quenching. 

Similar  to  photon-molecule  collisions,  electron-molecule  collisions  may  be  elastic  or 
inelastic.  With  inelastic  collisions,  the  molecule  is  excited  to  a  higher  energy  level  and  the 
electron  loses  kinetic  energy.  Photon  emissions  occur  when  the  molecule  drops  from  the 
excited  state  to  the  ground  state.  The  primary  interactions  considered  for  this  study  are  photon- 
molecule  inelastic  scattering  and  flourescense  scattering. 
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Figure  II-5  Molecular  processes  induced  by  incident  radiation. 
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III.     DAYGLOW  EMISSION  FEATURES 
A.    Molecular  Nitrogen 

Figure  IQ-1  is  an  N2  energy  level  diagram  including  the  two  transitional  bands  of  interest 
in  the  1800-3400  A  range.  The  states  are  indicated  by  molecular  notation,  and  the  transitions 
are  identified  by  name.  Since  these  are  forbidden  transitions,  the  excitation  is  from  electron- 
molecule  inelastic  scattering. 

1.  Vegard-Kaplan  Bands 

The  Vegard-Kaplan  system  is  so  designated  because  Vegard  discovered  the  system 
in  1932,  and  Kaplan  determined  the  specific  wavelengths  involved  a  few  years  later.  The 
emissions  range  from  1250  A  to  5325  A  and  represent  the  A3ZU+->X,Z  +  transition.    In  this 

transition,  the  total  transformation  is  +->+,  the  electronic  inversion  is  u— >g,  AZ=  -1,  and  AA=0. 
From  the  molecule  selection  rules,  this  is  a  forbidden  transition,  yet  it  has  been  observed  as  an 
electric  dipole  emission.  A  thorough  discussion  of  the  band  system  is  given  in  Meier  (1990). 

A-state  nitrogen  molecules  are  produced  by  cascading  downward  from  the  B  and  C 
states  or  as  a  result  of  photoelectron  collisions.  Because  the  relationship  of  emission  bands  to 
photoelectron  flux  is  complicated  by  cascading  N2  molecules  and  quenching  effects,  the  V-K 
band  is  not  generally  used  to  infer  photoelectron  flux. 

2.  Lyman-Birge-Hopfield 

The  Lyman-Birge-Hopfield  system  is  the  name  given  to  the  A1ng-»X1Zg+  transition 
of  N2.  The  emission  bands  arising  from  this  transition  extend  from  1 160  A  to  3020  A  and  are 
due  to  both  electric  quadrupole  and  magnetic  dipole  transitions.    Unlike  the  V-K  system,  the 
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upper  state  of  the  LBH  system  can  only  be  populated  by  direct  photoelectron  collision 
excitation.  This  system  can,  therefore,  be  directly  used  to  infer  the  photoelectron  flux.  The 
LBH  bands  found  in  the  MUSTANG  data  are  observed  in  the  wavelength  region  from  1 800  A 
to  2400  A.  An  intensive  investigation  of  the  LBH  emissions  observed  during  the  1 990  flight 
was  performed  by  Mack  (1991). 

B.    Nitric  Oxide 

Nitric  oxide  (NO)  emissions  dominate  much  of  the  dayglow  data.  There  are  four  primary 
band  systems  in  the  range  1800  A-3400  A:  y,  5,  P,  and  e.  These  NO  transitions  are  generated 
by  photoexcitation.  The  P  bands  are  several  orders  of  magnitude  less  intense  than  the  other 
emissions,  so  they  have  been  neglected  in  this  analysis.  Figure  ni-2  is  an  NO  energy  level 
diagram  including  the  three  transitional  bands  of  interest  in  the  1 800  A-3400  A  range. 

1.  Delta  Bands 

The  C2n->X2TI  transition  band  system,  designated  5,  ranges  from  1600  A  to  4300  A 
(Barth,  1965).  Using  standard  notation,  the  upper  electronic  state's  vibrational  level  is 
designated  v1,  and  the  lower  state  vibrational  level,  v".  The  5  transition  includes  vibrational 
levels  corresponding  to  v-0,1,2,3,4  and  v'-0,l,...23.  Nitnc  oxide  will  dissociate  at  v*  >  0  for 
C2n,  so  v'  is  limited  to  v-  0  only. 

2.  Gamma  Bands 

The  y  bands  of  nitric  oxide  extend  from  1 850  A  to  61 50  A.  The  emissions  arise  from 
the  A2I+-»X2n  transition,  and  the  ground  state  has  v"=  0,...23  as  previously  stated.  The 
dissociation  energy  for  nitric  oxide  is  lower  than  the  energy  of  the  v  -  4  level  in  the  A  state,  so 
only  transitions  arising  from  v-  0,1,2,3  contribute  to  the  synthetic  y  spectrum. 
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3.     Epsilon  Bands 

The  nitric  oxide  e  system  is  the  result  of  molecules  transitioning  from  D2Z+— >X2n. 
The  wavelengths  of  emissions  range  from  1 800  A  to  2400  A.  D-state  molecules  populate 
vibrational  levels  v*=  0,1,2. 
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IV.     THE  EXPERIMENT 
A.    Introduction 

The  middle  ultraviolet  spectra  analyzed  in  this  thesis  were  obtained  from  the  NPS 
MUSTANG  spectrograph  flown  on  a  NASA  sounding  rocket  launched  on  March  19,  1992, 
from  White  Sands  Missile  Range,  New  Mexico.  The  MUSTANG  instrument  was  launched  on 
a  Terrier  boosted  Black  Brant  vehicle,  which  carried  the  payload  to  an  apogee  altitude  of 
320km.  The  spectrograph  field-of-view  was  pointed  in  the  anti-solar  direction  with  an 
observation-zenith- angle  of  1 00°  on  the  up-leg  and  90°  on  the  down-leg. 


B.     Instrument  Description 

The  MUSTANG  (Middle  Ultraviolet  Spectrograph)  shown  in  Figure  IV- 1,  covers  a 
wavelength  range  from  1850A  to  3400 A.  It  is  a  modified  1/8  m  Ebert-Fastie  spectrograph 

Speclrograph   Entrance  Slii 


l/8th  m  Telescope   Mirror 


llamamatsu  Image  Sensor 

Figure  IV- 1  Schematic  diagram  of  the  MUSTANG  instrument. 
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with  a  1/8  m  telescope  constructed  by  Research  Support  Instruments,  Inc.  The  entrance  slit  is 
5  mm  by  140  urn  which  yields  a  field-of-view  of  2.3°  by  0.07°.  A  10  A  wavelength  resolution 
is  provided  by  a  HYPERFTNE  (Inc.)  1200  lines/mm  grating  ruling.  The  detector  assembly 
consists  Of  an  ultraviolet-to-visible  ITT  image  intensifier  and  a  Hamamatsu  linear  array 
photodiode  detector.  The  image  intensifier  uses  a  quartz  input  window,  a  CsTe  photocathode, 
and  a  fiber  optic  output  window.  The  detector  is  an  array  of  5 1 2  photodiodes  of  area  5  mm  by 
50  urn  giving  analog  voltage  outputs  proportional  to  intensity  using  a  fixed  50  msec  integration 
period.  This  voltage  output  is  fed  serially  to  an  A/D  converter  and  then  to  rocket  telemetry  via 
a  FIFO  shift  register. 

Pre-launch  calibration/sensitivity  was  obtained  using  standard  techniques  with  a 
Deuterium  lamp  for  the  wavelength  range  1 800  A-2800  A  and  a  FEL  tungsten  lamp  for  the 
range  2400  A-3400  A.  The  calibration  curves  obtained  from  the  two  lamps  agree  to  within 
10%.  The  MUSTANG  instrument  experienced  no  significant  degradation  in  sensitivity  (less 
than  5%)  during  the  launch  sequence  based  on  post-flight  calibration. 

C.    Wavelength  Calibration 

The  wavelength  calibration  was  determined  using  the  emission  spectra  of  Mercury  and 
Platinum.  The  platinum  lamp  provides  25  spectral  emissions  in  the  wavelength  region  of 
interest.  The  design  of  the  MUSTANG  instrument  is  optimized  for  a  bandwidth  from  1850  A 
to  3400  A  and  a  dispersion  of  3. 13  A/pixel. 

The  pixel  positions  on  which  the  selected  platinum  emissions  fell  were  determined  by  a 
rough  approximation  (  ±  0.5  A)  of  the  centroid  of  each  spectral  line.  The  wavelengths  of  the 
25   emissions   were  plotted  versus  pixel  position.      Previous   works   (Anderson,    1990; 
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Chase,  1992)  used  a  linear  best  fit  calibration  procedure.     It  was  decided  to  attempt  a 
polynomial  best  fit  for  this  analysis.  The  polynomial  fit  was  determined  to  be: 

A=[l848.270  +  3.188Ar-(1.374X10-4)AT2]  A; 
where  TV  is  the  pixel  number. 

D.    Data  Collection 

Over  8000  spectra  were  collected  by  the  MUSTANG  instrument  during  the  rocket 
experiment.  These  spectra  were  averaged  into  65  altitude  bins,  disregarding  a  few  spectra 
which  were  contaminated  by  noise  spikes  from  the  instrument's  electronics.  Each  of  the  65 
averaged  spectra  covers  a  5  km  layer  that  is  labeled  by  its  central  altitude.  The  adjusted 
spectra  were  then  divided  by  instrument  sensitivity.  Each  spectrum  contained  512  intensity 
values.  The  first  two  values  were  set  to  zero,  because  the  uncertainty  in  the  instrument 
sensitivity  at  the  lower  wavelength  endpoint  is  extremely  large. 
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V.       DATA  ANALYSIS 
A.    Algorithm  Techniques 

1.  Synthetic  Spectra 

Computer  synthetic  spectra  written  by  Geary  (1986)  were  used  for  the  generation  of 
Nitric  Oxide  5,  e,  and  y  and  the  Nitrogen  Lyman-Birge-Hopfield  band  emissions.  The 
Nitrogen  Vegard-Kaplan  synthetic  model  was  written  by  Siskind  and  Barth  (1987).  All  of  the 
computer  code  used  in  this  study  were  written  in  Interactive  Data  Language  (IDL®),  Version 
3.1,  designed  by  Research  Systems,  Inc.  The  inputs  for  the  synthetic  spectrum  algorithms 
were  temperature  and  the  calibrated  wavelength  range  of  the  MUSTANG  instrument.  The 
output  is  an  array  of  512  intensities  corresponding  to  the  instruments  wavelength  range. 
Before  comparing  the  synthetic  spectra  with  the  data,  the  synthetic  spectra  were  convolved 
with  MUSTANG'S  instrument  response  function.  The  MUSTANG  instrument  response  was 
estimated  to  be  an  1 1  -element  symmetric  slit  function  (Clayton,  1 990;  Anderson,  1 990;  Mack, 
1991). 

2.  Fitting  Method 

The  technique  used  to  fit  the  synthetic  spectra  to  the  data  was  an  iterative  process  of 
modifying  a  linear  combination  of  each  band  emission  system  and  adjusting  separate  synthetic 
and  quantum  parameters  of  each  emission  as  necessary.  These  adjustable  parameters  are  the: 
1)  Intensity  Scale  Factor;  2)  Molecular  Oscillator  Strength;  and  3)  the  Franck-Condon  Factor. 

The  intensity  scale  factor  is  a  relative  weight  assigned  to  each  band  system.  It  is  a 
best  "initial  guess"  taking  into  account  the  combined  contributions  of  other  surrounding 
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emissions.  For  example,  the  entire  N2  LBH  model  was  assigned  a  specific  scale  factor  and 
then  summed  with  all  other  contributing  spectra. 

The  molecular  oscillator  strength  is  a  quantum  mechanical  component  for  the 
synthetic  models  which  quantifies  the  probability  of  a  molecule  absorbing  a  photon  in  the 
electronic  ground  state  and  transitioning  to  a  specific  vibrational  level  (v')  in  the  excited 
electronic  state.  This  mechanism  of  molecular  photon  absorption  is  discussed  in  the 
background  section.  For  instance,  in  the  N2  VK  synthetic  band  system,  at  a  temperature  of 
1 125  K,  all  v-  3  bands  were  assigned  a  scale  factor  of  .7  in  order  to  decrease  the  intensity  by 
30%  for  a  best  fit. 

The  Franck-Condon  factor  is  a  quantum  mechanical  element  which  quantifies  the 
probability  of  a  synthetic  emission  occurring  from  a  specific  vibrational  level  (v')  in  the  excited 
electronic  state  to  a  specific  vibrational  level  (v")  in  a  lower  or  ground  electronic  state.  For 
example,  within  the  NO  y  synthetic  band  system  the  y  (2,3)  emission  was  observed  to  be 
underestimated  by  50%  but  the  y  (2,6)  was  overestimated  by  35%  and  both  were  assigned  a 
specific  scale  factor  to  correct  the  discrepency. 

The  adjustable  synthetic  and  quantum  parameters  designated  as  the  intensity  scale 
factor,  molecular  oscillator  strength,  and  the  Franck-Condon  factor  will  be  referred  to  as  the 
"SQ  parameters"  for  the  rest  of  this  work. 

B.    Intensity  Scale  Factor  Determination 

Intensity  scale  factors  were  determined  by  visual  inspection  for  "goodness  of  fit".  An 
analysis  strategy  was  developed  to  take  advantage  of  various  attributes  of  NO  and  N2 
molecular  spectra.  For  example,  the  N2  LBH  and  NO  e  band  systems  are  more  prominent 
from  1900  A  to  2200  A  while  the  N2  VK  and  NO  y  band  systems  are  more  dominant  in  the 
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2200  A  to  3100  A  range.  A  time  intensive  iterative  approach  was  conducted  to  identify  each 
synthetic  band  system  which  had  a  discrepency  with  the  data.  An  initial  attempt  was  made  to 
assign  one  scale  factor  per  band  system.  However,  it  was  necessary  to  modify  individual 
emissions.  The  SQ  parameters  of  these  emissions  were  modified  to  achieve  a  best  fit  synthetic 
spectrum  to  the  rocket  data. 

C.  Nitric  Oxide  Self-Absorption 

A  self  absorption  model  using  Holstein  transmission  functions,  developed  by  Eparvier  and 
Barth  (1992),  was  applied  to  the  NO  (0,0),  (1,0),  (2,0)  y  resonance  band  emissions.  The 
generated  synthetic  models  were  then  added  to  other  contributing  spectra.  If  the  intensities  had 
a  discrepency  to  the  data,  their  respective  SQ  parameters  were  modified  for  a  best  fit. 

D.  NO  Column  Density  Determination 

The  NO  column  density  profile  determination  was  performed  in  tandem  with  the 
temperature  profile  and  NO  intensity  scale  factor  estimations.  The  NO  column  density  (poi) 
for  each  altitude  bin  was  calculated  from  the  following  relationship: 


Where:  sr  =  3.13 


'  A  N 
^  pixel  > 


rj  =  NO  intensity  scale  factor 


R 


u 


R  =  1X10 


f  photons  ^  [atom  *  pixel  *  sec^ 


Vcm  *sec.' 


V         photon        J 

The  k  term  is  the  MUSTANG  instrument  pixel  spacing.    The  dependence  on  instrument 
sensitivity  has  been  removed  from  the  rocket  data  spectra  to  yield  absolute  spectral  intensity 
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(r\)    in    units    of    Rayleighs    per    Angstrom     \~r\-        A    Rayleigh    (R)    is    defined 


as     1X10 


f  photons  ' 
vcm2  *  secy 


.    And  the  C,  term  is  a  conversion  factor  representing  one  photon 


emitted  from  each  NO  atom  and  detected  by  an  individual  pixel  of  the  MUSTANG  instrument 
per  unit  time. 

£.    Temperature  Estimations 

The  temperature  of  each  altitude  bin  was  estimated  by  isolating  an  observed  molecular 
band  relatively  free  of  other  contaminating  spectra.  The  standard  emission  used  below  1 60  km 
was  the  NO  y  (0,1)  at  2365  A.  Above  160  km  the  N2  VK  system  begins  to  contaminate  the 
NO  y  (0,1).  The  N2  VK  (0,5)  at  2605  A  was  used  for  estimating  a  temperature  profile  for  all 
of  the  upleg  and  downleg  altitude  bins  above  160  km. 
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VI.     RESULTS 
A.    Comparison  of  Data  and  Synthetic  Spectra 

The  averaged  spectra  examined  here  are  for  the  altitude  data  bins  1 1 0,  1 60,  and  1 95  km. 
The  primary  emission  features  observed  in  this  MUV  spectral  region  of  1 900  A  to  3 1 00  A 
include  the  N2  Lyman-Birge-Hopfield  (LBH),  N2  Vegard-Kaplan  (VK),  and  the  NO  Gamma 

(y),  Delta  (6),  and  Epsilon  (e)  bands.  All  known  emissions  in  this  range  are  used  in  the  fitting 
routines  except  OI 2972  A,  Oil  2470  A,  Om  2853  A,  and  NH  2143  A  atomic  lines. 

A  comparison  of  the  observed  spectra  with  the  synthetic  spectrum  model  in  the  spectrum 
for  all  three  of  the  altitude  ranges  is  shown  in  Figures  A-l  to  A-22  in  Appendix  A.  In  order  to 
illustrate  the  basic  technique  and  results  of  the  fitting  technique  employed  for  this  work,  the 
160  km  data  set  was  used  as  a  sample  (Figures  A-l  to  A- 10).  The  160  km  data  set  was  chosen 
because  the  intensities  of  the  NO  and  N2  band  emission  systems  were  approximately  equal. 

Figure  A-l  illustrates  the  160  km  NO  e  (v-0)  emissions  for  the  range  1900  A  to  2130  A. 
The  N2  LBH  (v-3)  emissions  in  the  same  wavelength  range  are  shown  in  Figure  A-2.  Both  of 
these  synthetic  spectra  are  given  to  show  the  individual  contributing  emissions  of  the  five  band 
systems.  To  show  a  complete  band  system,  the  sum  of  all  the  N2  LBH  emissions  is  shown  in 
Figure  A-3.  Figure  A-4  presents  the  summation  of  the  N2  LBH  and  NO  e  band  systems. 
Continuing  with  the  step  by  step  additions  of  the  other  three  band  systems,  the  result  is  shown 
in  Figure  A-5. 

The  individual  emissions  were  used  to  identified  discrepancies  between  the  data  and 
synthetic  spectra.  This  was  accomplished  by  overplotting  several  emissions  in  the  wavelength 
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range  in  question.  Once  an  element  was  identified  as  the  erroneous  contributor  (e.g.,  exceeded 
or  underestimated  observed  intensities),  its  respective  SQ  parameters  were  modified  then 
resummed  for  a  best  fit. 

Figure  A-5  reveals  several  interesting  features  of  the  1 60km  data  set.  Note  that  only 
minimum  intensity  scale  factors  are  applied.  One  can  visually  identify  the  synthetic  spectrum 
exceeding  the  observed  data  at  approximately  1910  A,  1980  A,  2050  A,  and  2095  A.  The 
synthetic  spectrum  can  also  be  seen  lagging  at  2030  A  and  2025  A. 

Figure  A-6  presents  the  data  versus  wavelength  compared  with  the  theoretical,  but 
modified  oscillator  strengths  to  v'  emissions  have  also  been  applied.  Notice  that  the  model 
underestimated  the  observed  intensity  at  1910  A  and  1930  A.  To  correct  this,  Figure  A-7 
shows  how  modifying  the  Frank-Condon  factors  for  each  erroneous  (V,  v")  emission  adjusted 
the  theoretical  fit  closer  to  the  real  data  Figures  A-8,  A-9,  and  A- 10  repeats  this  process  and 
displays  the  entire  160  km  data  set  for  1900  A  to  3100  A.  In  addition,  all  of  the  altitude  data 
bins  at  1 10,  160,  and  195  km  are  illustrated  in  Figures  A-l  1  to  A-22. 

The  resulting  quantum  molecular  oscillator  strengths  (v1)  and  Franck-Condon  factors  (v") 
to  achieve  a  best  fit  for  each  band  emission  for  all  three  altitude  data  sets  are  listed  in  Tables  1 
through  3  in  the  Appendix  section. 

B.    Nitric  Oxide  Bands 

The  nitric  oxide  synthetic  spectra  were  substantially  contaminated  by  molecular  nitrogen  at 
160  km.  The  y  (0,1)  is  an  intensity  standard  at  lower  altitudes  but  is  interfered  with  by 
surrounding  N2  Vegard  Kaplan  emissions. 

A  successful  attempt  was  made  to  standardize  the  NO  Gamma  (y),  Delta  (5),  and  Epsilon 
(e)  bands  by  weighing  all  three  systems  by  the  same  intensity  scale  factor.     This  was 
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accomplished  by  modifying  key  values  of  the  SQ  parameters  for  each  band  emission.  The  NO 
Delta  (5),  and  Epsilon  (e)  bands  were  the  most  difficult  to  model  due  to  contaminating  N2 

Lyman-Birge-Hopfield  (LBH)  emissions  in  the  1900  A  to  2120  A  range.  The  NO  Gamma  (y) 
band  system  was  also  contaminated  with  N2  VK  emissions  in  the  2200  A  to  3 1 00  A  range,  but 

to  a  lesser  degree. 

C.  Molecular  Nitrogen  Bands 

At  160  km  the  Vegard  Kaplan  bands  begin  to  dominate  over  all  other  emissions  with  the 
exception  of  the  y  (1,0).  Initial  analysis  indicated  two  discrepencies  of  the  synthetic  spectra 
concerning  these  bands.  First,  a  negative  3.13  A  shift  (approximately  1  pixel  resolution)  of  the 
VK  synthetic  versus  data  spectra  To  correct  this  discrepency,  a  modification  to  the  VK 
synthetic  algorithm  (VKTERM.PRO)  of  the  initial  molecular  constants  was  performed.  The 
TO  molecular  constant  was  change  from  50203.6  to  50169.0  and  the  WXEE  molecular 
constant  was  changed  from  14.32  to  11.45.  And  the  second  discrepency  was  the  observed 
relative  intensities  of  the  VK  (0,5),  (0,6),  (0,7)  bands.  The  synthetic  predicted  a  constant  ratio 
of  these  three  emissions  over  all  temperature  inputs.  However,  the  observed  data  showed 
varying  scale  factors  which  were  matched  by  modifying  the  synthetic  model  and  have  been 
recorded  in  Appendix  D,  Tables  B-l  through  B-3. 

D.  Nitric  Oxide  Self-Absorption 

A  self  absorption  model  was  made  of  the  theoretical  self  absorption  versus  non-self 
absorption  mechanisms  of  the  NO  (0,0),  (1,0),  (2,0)  y  resonance  band  emissions.  The  Holstein 
transmission  functions  applied  for  this  model  generated  mixed  results.  The  NO  (0,0)  y 
resonance  band  was  consistently  overestimated  at  all  three  altitudes  analyzed.   The  NO  (1,0), 
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(2,0),  y  resonance  bands  also  had  to  have  modifications  to  their  respective  SQ  parameters,  but 
the  results  did  not  indicate  any  consistency.  Figure  A- 11  through  A- 14  includes  the  NO  y 
resonance  bands  at  1 1 0  km. 

E.  Nitric  Oxide  Column  Density  Profile 

The  nitric  oxide  density  profile  inferred  from  the  data  is  determined  from  all  of  the 
recorded  altitude  bins  and  is  shown  in  Appendix  B,  Figure  B-l  and  Table  B-l.  These 
estimated  densities  are  consistent  with  other  observations  during  similar  solar  cycle  intensities. 

F.  Temperature  Profile 

A  temperature  profile  estimated  from  the  data  is  determined  from  all  of  the  recorded 
altitude  bins  and  is  shown  in  Appendix  C,  Figure  C-l  and  Table  C-l.  These  estimated 
temperatures  are  consistent  with  other  observations  during  similar  solar  cycle  intensities. 
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VII.    CONCLUSIONS 

A.  Summary  of  Findings 

In  this  work,  spectra  from  the  Earth's  ionospheric  dayglow  obtained  from  a  sounding 
rocket-launched  experiment  were  analyzed.  Computer  generated  spectra  were  created  for  the 
NO  y,  5,  and  s  bands  and  the  N2  LBH  and  VK  bands.   These  synthetic  spectra  were  fitted  to 

the  observed  dayglow  data.  An  overall  synthetic  fit  agreed  well  with  the  observed  spectra. 
Intensity  scaling  factors,  molecular  oscillator  strengths  and  Franck-Condon  factors  were  used 
and  modified  to  obtain  a  best  fit.  Synthetic  emission  bands  that  were  observed  to  be  in 
discrepency  with  the  observed  data  were  visually  identified  and  corrected  to  obtain  a  best 
possible  fit. 

The  NO  intensity  scaling  factors  were  used  to  infer  a  NO  column  density  profile.  A 
temperature  profile  was  estimated  from  the  synthetic  spectra  which  best  fit  the  data.  The 
synthetic  fit  also  confirmed  the  theoretical  self  absorption  versus  non-self  absorption  processes 
of  the  NO  (0,0),  (1,0),  (2,0)  y  resonance  band  emissions.  NO  self  absorption  is  a  necessary 
consideration  of  any  atmospheric  nitric  oxide  analysis  scheme. 

B.  Recommendations  For  Further  Research 

The  NPS  MUSTANG  will  be  part  of  an  operational  remote  sensing  spacecraft.  It  will  be 
launched  on  the  Advanced  Research  and  Global  Observation  Satellite  (ARGOS)  in  late  1995. 
Considering  the  time  intensive  nature  of  a  visual  analysis  by  comparing  the  observed  data  with 
synthetic  spectra,  a  high  speed  central  processor  is  strongly  recommended.  Once  an  intensive 
catalog  of  synthetic  emissions  (which  are  in  discrepency  with  observed  airglow)  is  complete, 
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appropriate  modifications  to  the  fitting  programs  would  add  to  the  reliability  of  MUSTANG  as 
a  surveyor  of  ionosperic  parameters  for  national  defense  agencies. 
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VIII.  APPENDIX  A 

Appendix  A  contains  the  figures  of  rocket  data  versus  wavelength  and  comparing  the 
data  with  the  generated  theoretical  spectrums  as  referred  to  in  the  results  section. 
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IX.     APPENDIX  B 

Appendix  B  contains  Figure  B-l  which  shows  the  NO  column  density  profile  versus 
altitude  and  Table  B-l  which  lists  each  data  point  of  the  NO  column  density  profile  as  referred  to  in 
the  results  section. 
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X.       APPENDIX  C 

Appendix  C  contains  Figure  C-l  which  shows  the  temperature  profile  versus  altitude  and 
Table  C-l  which  lists  each  data  point  of  the  temperature  profile  as  referred  to  in  the  results  section. 
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XL     APPENDIX  D 

Appendix  D  contains  the  listings  for  minimum  scale  factors,  modified  oscillator  strengths, 
and  modified  Franck-Condon  factors  (Best  Fits)  corresponding  to  the  data  bins  110,  160,  and 
195km.  and  are  shown  in  Tables  D-l  thru  D-4. 

Table  D-l  lists  the  intensity  scale  factors  for  the  N2  LBH  and  VK  synthetic  spectra,  and 

the  NO  y,  8,  and  e  synthetic  spectra. 

Tables  D-2,  D-3,  and  D-4  contain  the  fractions  necessary  to  modify  each  band  emission  for 
a  best  fit.  If  a  fraction  is  entered  for  a  v'  row,  it  is  equivalent  to  applying  that  percentage  to  the 
molecular  oscillator  strength.  If  a  fraction  is  entered  for  a  v"  row,  it  is  equivalent  to  applying  that 
percentage  to  the  molecular  oscillator  strength  and  the  Franck-Condon  factor.  If  any  entry  is 
blank,  assume  the  factor  is  1.0. 
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Table  D-l 


Synthetic  Scale  Factors 


Altitude  (km): 

110 

160 

195 

Temperature  (Kelvin): 

200 

650 

1125 

VK 

7.59E+03 

1.46E+04 

1.59E+04 

LBH 

4.46E+04 

2.50E+04 

1.46E+04 

NO 

3.51E+09 

5.46E+08 

2.79E+08 

NO  Column  Density: 

9.08E+15 

1.22E+15 

5.13E+14 

Synthetic  Self  Absorption       Scale  Factors 

0,0 

1.00 

0.96 

0.75 

1,0 

0.82 

0.97 

0.72 

2,0 

0.62 

0.48 

0.65 
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Table  D-2 

Altitude:     110  km 

Temp: 

200  K 

v' 

v" 

VK 

LBH 

GAM 

EPS 

DEL 

V 

v" 

VK 

LBH 

GAM 

EPS 

DEL 

0 

0.98 

0.21 

4 

0 

0.12 

5 

0.45 

1 

0.58 

0.12 

7 

0.32 

2 

0.54 

8 

0.03 

3 

0.86 

0.79 

0.62 

10 

0.32 

4 

0.66 

11 

0.2 

5 

0.77 

0.95 

12 

0.3 

6 

13 

0.3 

7 

14 

0.3 

1 

0.66 

15 

0 

0.83 

5 

1 

5 

2 

6 

3 

0.45 

0.6 

7 

4 

0.74 

0.6 

8 

5 

0.68 

9 

0.27 

6 

0.66 

10 

0.01 

7 

0.55 

11 

8 

0.86 

12 

0.3 

2 

0.66 

13 

0.3 

0 

14 

1 

1 

0.74 

15 

2 

6 

0.23 

3 

0.41 

8 

0.59 

4 

0.55 

9 

0.05 

5 

0.14 

0.66 

7 

6 

0.55 

6 

0.45 

7 

0.59 

10 

0.27 

8 

0.55 

11 

0.01 

9 

0.03 

8 

10 

0.03 

9 

3 

0.5 

7 

0.68 

0 

0.58 

9 

0.01 

1 

11 

0.01 

2 

12 

0.14 

3 

10 

4 

13 

0.05 

5 

0.45 

n 

6 

0.09 

9 

7 

12 

0.01 

8 

0.45 

14 

0.27 

9 

13 

10 

.04 

9 

0.45 

11 

0.04 

11 

0.23 

12     | 

0.04 

13 

0.45 
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Table  D-3 

Altitude:     160  km 

Temp: 

350  K 

v' 

v" 

VK 

LBH 

GAM 

EPS 

DEL 

v' 

v" 

VK 

LBH 

GAM 

EPS 

DEL 

0 

0.54 

4 

0 

0.81 

0.38 

4 

0.74 

1 

0.72 

0.86 

0.28 

5 

0.39 

2 

0.46 

0.38 

7 

0.27 

3 

0.66 

0.77 

0.56 

8 

0.01 

4 

0.88 

0.21 

10 

0.27 

5 

0.81 

11 

0.42 

6 

1 

0 

12 

0.86 

7 

0.63 

13 

0.6 

1 

14 

0.4 

0 

0.83 

1 

0.77 

5 

2 

6 

0.39 

3 

0.39 

0.71 

0.58 

9 

0.23 

0.56 

4 

0.98 

0.78 

0.42 

10 

0.01 

5 

0.62 

0.79 

0.42 

12 

0.62 

6 

0.75 

13 

7 

0.51 

0.79 

6 

0.02 

8 

0.56 

0.46 

8 

0.51 

2 

0.1 

9 

0.04 

0 

7 

1 

6 

0.39 

2 

0.58 

7 

0.55 

3 

0.35 

0.52 

10 

0.23 

4 

0.31 

11 

0.01 

5 

0.01 

0.67 

8 

6 

0.39 

7 

0.01 

7 

0.35 

9 

8 

0.23 

7 

0.7 

9 

0.59 

9 

0.01 

10 

0.46 

11 

0.01 

3 

12 

0.12 

0 

10 

1 

0.5 

13 

0.04 

2 

11 

3 

8 

0.01 

4 

10 

0.01 

5 

0.35 

12 

0.04 

6 

0.08 

12 

7 

9 

8 

0.39 

12 

0.01 

9 

14 

0.23 

10 

0.79 

13 

11 

0.7 

9 

0.39 

12 

0.14 

11 

0.2 

13 

13 

0.39 

59 


Table  D-4 

Altitude:     195  km 

Temp: 

1125  K 

v' 

v" 

VK 

LBHl 

GAM 

EPS 

DEL 

v' 

v" 

VK 

LBH 

GAM)    GPS 

DEL 

0 

0.53 

0.19 

4 

0.01 

0 

0.71 

0.1 

4 

1 

0.59 

0.48 

0.01 

5 

2 

0.61 

0.06 

0.14 

7 

3 

0.77 

0.42 

0.42 

0.1 

8 

4 

0 

0.71 

0.42 

10 

5 

0.77 

11 

0.68 

6 

12 

0.8 

7 

0.88 

13 

0.69 

1 

0.52 

14 

0.05 

0 

0.65 

15 

1 

0.49 

5 

0.33 

2 

5 

0.01 

3 

0.42 

0.06 

6 

0.09 

4 

0.91 

0.58 

0.06 

7 

5 

0.61 

0.58 

8 

6 

0.1 

9 

0.01 

7 

0.53 

10 

0.01 

8 

0.58 

11 

2 

0.42 

0.52 

12 

0.76 

0 

0.48 

13 

0.11 

1 

6 

0.28 

0.47 

2 

0.45 

9 

0.01 

3 

0.51 

7 

0.19 

4 

6 

1.03 

5 

0.01 

8 

0.33 

6 

0.47 

0.01 

10 

0.19 

7 

0.47 

11 

0.01 

8 

0.79 

8 

0.28 

9 

9 

0.05 

10 

0.8 

9 

0.01 

3 

0.65 

10 

0.01 

0 

8 

0.47 

1 

11 

0.14 

2 

8 

3 

10 

0.01 

4 

12 

0.01 

5 

0.01 

12 

0.01 

6 

0.19 

9 

7 

12 

8 

0.09 

14 

9 

0.33 

13 

0.19 

10 

0.67 

9 

11 

0.95 

11 

0.01 

12 

0.58 

13 
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